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Abstract 

We study the dynamics of island nucleation in the presence of ad- 
sorbates using kinetic Monte Carlo simulations of a two-species growth 
model. Adatoms (A-atoms) and impurities (B-atoms) are codeposited, 
diffuse and aggregate subject to attractive AA- and AB-interactions. 
Activated exchange of adatoms with impurities is identified as the key 
process to maintain decoration of island edges by impurities during 
growth. While the presence of impurities strongly increases the island 
density, a change in the scaling of island density with fiux, predicted 
by a rate equation theory for attachment-limited growth [D. Kandel, 
Phys. Rev. Lett. 78, 499 (1997)], is not observed. We argue that, 
within the present model, even completely covered island edges do not 
provide efficient barriers to attachment. 
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The effect of impurities on crystal growth has been a long-standing con- 
cern in surface science Recent interest in this subject has been fueled by 
the prospect of using adsorbates as "surfactants" to improve the quality of 
epitaxially grown films [0. A key mechanism which dramatically increases 
the (detrimental or beneficial) effects of adsorbates is their tendency to dec- 
orate the step edges by preferentially attaching there P, ^. If this induces 
an additional energy barrier which an adatom has to overcome to be incor- 
porated at the step edge, the scaling exponent x iii the relation 

N ~ {F/DY (1) 

between the island density A^, the deposition flux F and the adatom diffusion 
coefficient D may be affected. Assuming a critical (= largest unstable) island 
size i* , standard rate equation theory yields in two dimensions the expression 
ii 

i* 

^ ~ i* + 2 

while strong incorporation barriers imply 

Thus inspection of the flux dependence of the island density could reveal 
whether or not the island edges are efficiently passivated by the adsorbates 
provided the critical island size is only slightly changed by the presence of 
adsorbates. Markov |§] has pointed out that the change from (0) to @ can 
be understood as a transition from a diffusion limited to a kinetically limited 
growth regime. 

In this paper we report on an extensive numerical study of a two-species 
growth model aimed at answering the following two questions: First, what 
kinetic processes are necessary to generate an island morphology in which 
the impurities decorate the island edges? Second, how does the presence of 
impurities affect the island density, and what role does the passivation of 
island edges play in this context? 

We employ a full diffusion solid-on-solid model of epitaxial growth with 
two surface species, which we denote A and B. We suppose that particles of 
type A correspond to the growing material and particles of type B represent 
the impurities. Simulations start on a fiat substrate composed only of A 
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atoms. Both species are randomly deposited with generally different fluxes 
Fa and Fb- Here, we restrict ourselves to codeposition with Fa = Fb- 
Similar results for predeposition of adsorbates (B atoms) as well as a study 
of the effect of the concentration of adsorbates will be presented elsewhere. 
The migration of each surface atom is modeled as a nearest-neighbor hopping 
process with rate Ru = kc exp(— E'/j/Zc^T), where = 10^'^ Hz is an adatom 
vibration frequency, Ed is the hopping barrier, T is the substrate temperature 
and kB is Boltzmann's constant. The hopping barrier is the sum of a term 
from the substrate E'sub and a contribution from each lateral nearest neighbor 
Et^. Contributions depend on local composition: For each term we have the 
four possibilities AA, AB, BA and BB. Thus, the hopping barrier of adatom 
X {A or B) is 

E^n- {nlE^Z + nrE^^)^ (4) 

Y=A,B 

where E^ is the hopping barrier for a free X adatom on a substrate atom y , 
nl^ is equal to one if the substrate atom is of type Y and zero otherwise, n^^ 
is the number of nearest-neighbor X-Y pairs, and E^^ is the corresponding 
contribution to the barrier (symmetric in X and Y). Lateral interactions 
between impurity atoms are neglected {E^^ = 0). 

In the simulations reported here we used E^ — 0.8 eV, E^ = 0.1 
eV, E^J^ = 1.0 eV, E^^ = 0.1 eV, and the substrate temperature T = 
500 K. Other parameters for lateral interactions were varied. The system 
sizes rang ed from 300x300 to 500x500. The low values of E^^^ and Eg^ 
ensure that atoms deposited on top of an impurity immediately descend to 
the substrate. We use fluxes in the interval from 0.00025 ML/s to 0.25 ML/s. 
In the case of homoepitaxy, different values of E^^ correspond to different 
sizes of critical nucleus i*. There are two different situations E^^ > E^^ 
and E^^ < E^^. In equilibrium at a low temperature, the former case leads 
to the formation of islands composed inside mainly of A atoms with B atoms 
bounded near the edges, whereas in the latter case it is energetically more 
favorable when B atoms are inside the island. We observed, however, that 
for our parameters, growth always leads to intermixing of A and B atoms, 
in both cases > E^^ and E^^ < E^^. 

Thus the energetic bias favoring segregation is not sufficient to obtain 
configurations with impurities mostly at island edges. To achieve this, we 
have to introduce an additional thermally activated process of exchange of 
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an A atom approaching an island edge covered with an impurity. This is a 
similar process to that introduced in Ref. ||^ for simulation of homoepitaxy 
on Si(OOl) with predeposited hydrogen. Here, we allow the exchange of an 
A atom with an impurity when it has before the exchange process at most 
one bond to another A atom in a nearest-neighbor position. The rate of 
the exchange process is set to fcex = kQexp^—Eex/ksT), where the activation 
barrier E^x is taken to be independent of whether or not the A-atom has a 
nearest neighbor bond. The exchange barrier E^x was varied from 0.8 eV (the 
diffusion barrier of free adatoms) to 2 eV. For low Ec^ the impurities float 
on the island edges, whereas for large Eex they tend to get trapped inside an 
island. In order to obtain decorated island edges, i?ex has to be lower than a 
certain value which in the present case is about 1.2 eV; in the following we 
set £'ex = 1 eV. 

In Fig. la and lb we show examples of typical configurations with the 
same partial coverage of both species 9a = Ob = 0.1 ML (i.e. total coverage 
e = Oa+Ob = 0.2 ML) obtained by codeposition with fluxes Fa = Fb = 0.004 
ML/s. They illustrate the effect of varying the relation between E^^ and 
E^^ for E^^ = 0.3 eV (a typical configuration for homoepitaxy in shown 
in Fig. Ic). For E^^ = 0.2 eV some gaps can be seen in the impurity layer 
surrounding the islands, while for E^^ = 0.4 eV the decoration is complete. 
If we assume that the decorated island edge is in equilibrium with an ideal 
impurity gas of coverage 9b, then a simple detailed balance argument shows 
that the concentration of uncovered edge sites is (1+6'Be^" /ksT^-i^ Thus for 
T = 500 K, 9b = 0.1 and E^^ = 0.2 eV about one tenth of the edge sites are 
uncovered, in accord with the visual inspection of Figure 1. Correspondingly, 
simulations with E^^ = 0.1 eV show that large portions of the step edges 
remain uncovered. 

It can be seen from Figure 1 that the island density increases with in- 
creasing E^^. A quantitative evaluation of the dependence of island density 
on flux, coverage and E^^ is given in Figure 2. For comparison we show 
also results for homoepitaxial growth without impurities. We would like to 
point out three noteworthy features of these data. First, increasing the in- 
teraction energy between adatoms and impurities dramatically increases the 
island density. Second, the increase is essentially independent of flux, i.e. 
the scaling exponent x iii Eq.(|l|) is hardly affected by the impurities; we find 
X = 0.54 for E^^ = 0.2, x = 0.45 for E^^ = 0.4, and x = 0.54 for homoepi- 
taxial growth. Third, for large fluxes and large values of E^^, the island 
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density continues to increase with coverage at least up to 6' = 0.2 ML. 

While a detailed discussion of these effects will be left to an extended pub- 
lication, here we provide a simple argument to explain why the scaling expo- 
nent (0) predicted for passivated islands is not observed in our simulations. 
The mechanism described by Kandel relies on the rate of attachment of 
an adatom to an island, S, to be much smaller than the diffusion rate D on 
the terrace. It is useful to introduce the probability p = S/{D + S)^S/D 
for an adatom to attach during an encounter with an island edge. In the 
kinetically limited growth regime, characterized by @, an adatom typically 
has to visit many islands before being captured. To see when this is the 
case, consider an adatom diffusing through an array of islands of density 
and linear size R ~ \J6a/N . During n diffusion steps the adatom explores 
a region of area ~ n, which contains ~ Nn islands and ~ NnR edge sites. 
The adatom is captured when NnRp ^1. It follows that the condition for 
encountering many islands prior to capture reads 

p<^R-^ ^[N/e~A- (5) 

For a given value of p this places a lower bound on the island density (and, 
according to (|l]), a lower bound on the flux) beyond which a crossover to 
conventional, diffusion limited scaling will occur. 

Let us estimate p for a completely decorated island. An adatom ap- 
proaching the island has two options: It can return to the terrace with rate 
/coe~^^sub'^'^" )/fc-BT^ Qj, exchange with the impurity (and thus join the 

island) with rate fcoc"'^™'''^^"^. The capture probability p is the probability 
that the second process occurs before the first, and is given by 

With the standard values E^^ = 1 eV, E^^^ = 0.8 eV and E^^ > 0.2 eV it 
can be seen that the condition (^ is never satisfied. To significantly decrease 
p, one would have to either increase the exchange barrier E^^ relative to the 
diffusion barrier E^, which leads to intermixing of impurities and adatoms 
within the island, or to decrease the AB- interaction strength E^^, which 
would increase the number of gaps in the impurity layer and invalidate our 
assumption of complete decoration. In fact the appearance of E^^ in the 
exponent of (^ implies that even completely decorated island edges are not 
efficiently passivated in our model, because the B-atoms covering the edge 
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still have "bonds" available with which to attract A-atoms and keep them 
near the edge long enough for an exchange to occur. 

In conclusion, we have shown that the attachment-limited growth regime 
characterized by the scaling exponent is not easily realizable in a two- 
species growth system with isotropic nearest neighbor interactions. It re- 
mains to be understood what microscopic mechanism, other than passivation 
of island edges, is responsible for the strong increase in island density shown 
in Figure 2. The most obvious possibility is that the impurities reduce the 
mobility of adatoms by temporarily trapping them, thus effectively reducing 
the diffusion coefficient D in (0). A detailed analysis of this mechanism will 
be presented elsewhere. 

Acknowledgements. This work was supported by Volkswagen- Stiftung 
and by the COST project P3.130. 



6 



References 

[1] N. Cabrera and D.A. Vermilyea, in "Growth and Perfection of Crystals", 
ed. by R. Doremus, B. Roberts and D. Turnbull (Wiley, New York 1958), 
p. 393. 

[2] D. Kandel and E. Kaxiras, The surfactant effect in semiconductor thin 
film growth (preprint , |cond-mat /990 1 1 77| ) . 

[3] H.A. van der Vegt, H.M. van Pinxteren, M. Lohmeier, E. Vlieg and 
J.M.C. Thornton, Phys. Rev. Lett. 68, 3335 (1992). 

[4] M. Kalff, G. Comsa and Th. Michely, Phys. Rev. Lett. 81, 1255 (1998). 

[5] S. Stoyanov and D. Kashiev, in "Current Topics in Materials Science", 
ed. by E. Kaldis (North Holland, Amsterdam 1981), p. 69. 

[6] J. A. Venables, G.D.T. Spiller and M. Hanbiicken, Rep. Prog. Phys. 47, 
399 (1984). 

[7] D. Kandel, Phys. Rev. Lett. 78, 499 (1997). 
[8] I. Markov, Phys. Rev. B 56, 12544 (1997). 

[9] P. Smilauer, K. Mizushima and D. D. Vvedensky, Phys. Rev. Lett. 81, 
5600 (1998). 



7 



Figure Captions 



Fig. 1: Examples of configurations for flux Fa — 0.004 ML/s, adatom 
interaction energy E^^ = 0.3 eV, and coverage 9a = 0.1 ML: (a) codeposi- 
tion with Fb = Fa and E^^ = 0.2 eV, (b) codeposition with Fb = Fa and 
E^^ = 0.4 eV, (c) homoepitaxy {Fb = 0). We show only 50 x 50 sections of 
a larger simulation box. 

Fig. 2: Averaged island density as function of fiux Fa for several values 
of the total coverage 9 = 9a + Ob and different energy barriers: E^^ = 0.2 
eV - open symbols, E^^ = 0.4 eV - filled symbols. The adatom interac- 
tion energy E^"^ — 0.3 eV and the exchange barrier Eex = 1 eV are fixed, 
and the impurity flux Fb — Fa- The behavior in the absence of impurities 
(homoepitaxy, Fb — 0) is shown for comparison. 
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Figure 1: 
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Figure 2: 



10 



